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1.  INTRODUCTION 


The  research  carried  out  under  Grant  N00014-8D-J-1109,  which  began  1  October 
1988  and  is  the  successor  to  Contract  N00014-S4-K-0574,.is  primarily  in  the  field  of 
nonlinear  acoustics.  The  broad  goal  is  to  determine  the  laws  of  behavior  of  finite- 
amplitude  sound  waves,  especially  to  find  generalizations  of  the  known  laws  of  lin¬ 
ear  acoustics.  This  report  is  the  second  annual  report  under  the  Grant  and  covers 
the  12-month  period  ending  30  September  1990.  The  previous  report  (first  annual 
report’)  is  referred  to  herein  as  89-6.* 

The  following  persons  participated  in  the  •;esearch: 

Graduate  students 

C.  E.  Bradley,  M.S.  student  in  Mechanical  Engineering 
J.  A.  Ten  Cate,  Ph.D.  student  in  Mechanical  Engineering 

Undergraduate  student 
M.  J.  Van  Doren,  senior  in  Mechanical  Engineering 


Senior  personnel 

M.  F.  Hamilton,^  Mechanical  Engineering  Department,  The  University  of  Texas 
at  Austin 

J.  Naze  Tjptta,^  Research  Fellow,  on  leave  from  Mathematics  Institute,  Univer¬ 
sity  of  Bergen,  Norway 

W.  M.  Wright,  Consultant,  Physics  Department,  Kalamazoo  College,  Michigan 
D.  T.  Blackstock,  principal  investigator 


"iNmnbcr.s  given  in  t  his  stylo  roler  lo  ileins  in  liio  CliroiKjiugiuii  iJiliiioginpiiy  gi\on  at  tiie  oinl  of 
this  ri'port,  o.g.,  89-0  iiionii.s  the  .si.'ctli  entry  in  tiie  li.st  for  1989. 

Uliiniilton  received  no  direct  .support  from  llie  Grant.  However,  he  i.s  co-snpervi.sor  of 'leii  Cato’s 
IMi.l).  re.soarch,  wliicli  i.s  de.srrihed  in  Project  IJ  below. 

UMtIiongli  Nax.e  Ijotta  received  no  direct  support  from  Hie  Giant,  she  jiroviiled  sni)ervi.->ion  for 
the  tlieoretical  portion  of 'I'en Cate’s  I’h.l).  re.searcli,  wliicli  is  di'scrilied  in  i’rojort  11  below. 
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2.  PROJECTS 


The  following  projects  were  active  during  the  report  period: 

A.  Propagation  in  a  Periodic  Waveguide 

B.  Scattering  of  Sound  by  Sound 

C.  Laboratory  Upgrade 

D.  Miscellaneous 

Primary  effort  was  devoted  to  Projects  A  and  B. 


2.1  Propagation  in  a  Periodic  Waveguide 


This  is  Bradley’s  project.  Reported  last  year  were  theoretical  and  experimental 
results  for  acoustic  Bloch  wave  propagation  in  a  rectangular  waveguide  loaded  with 
rigidly  terminated  side  branches  (see  Sec.  2.2  in  89-6).  Finite-amplitude  as  well  as 
small-signal  waves  were  investigated.  The  primary  finding  for  small-signal  waves  was 
that  the  previously  derived  dispersion  relation  (Eq.  2.1  in  89-6)  showed  very  good 
agreement  with  Bradley’s  measurements  of  attenuation  and  dispersion.  Predictions 
from  the  nonlinjar  theory,  which  was  based  on  Korpel’s  progressive  spectral  wave 
equation,^  showed  good  qualitaiive  agreement  with  some  early  measurements  of  the 
second  harmonic  behavior,  but  poor  quaniiiaiivt  agreement  (89-4).  One  of  the  prob¬ 
lems  was  that  the  theoretical  predictions  failed  to  explain  the  strong  second  harmonic 
signal  that  was  measured  very  close  to  the  source  (it  was  determined  that  the  source 
itself  was  not  to  blame).  A  desire  to  improve  the  thcorc'tical  model  led  to  this  year’s 
work,  which  began  with  a  reexamination  of  the  linear  theory.  The  result  was  a  better 
understanding  of  the  structure  of  linear  Bloch  waves  in  a  periodic  waveguide  and,  in 
turn,  the  development  of  a  more  successful  nonlinear  theory. 

The  review  of  the  linear  theory  proved  to  be  quite  illuminating.  The  detailed 
structure  of  Bloch  waves  was  worked  out  and  a  new  expiession  of  the  Bloch  wave 
function,  referred  to  here  as  the  convolution  <  q)re.ssion,  was  derived.  Bloch  wave 
functions  can  be  expressed  in  a  variet}'  of  wa>'.  '  but  the  convolution  expression  aids 
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qualitative  understanding  of  the  structure  of  the  Bloch  wave  and  is  ideally  suited  to 
the  nonlinear  analysis.  Let  h  be  the  (spatial)  period  of  the  structure,  and  let  t/>(2) 
represent  the  pressure  wave  function  in  a  single  cell  — /i/2  <  z  <  /i/2  of  the  structure. 
It  is  found  that  the  Bloch  wave  function  p{z)  for  the  complete  structure  is  expressed 
as  the  convolution  of  xl){z)  with  a  phase  weighted  lattice  function, 

+00 

p{z)  =  tP{z)*  Yi  (2.1) 

n=-“00 


where  q  is  the  wave  number  of  the  Bloch  wave.  The  convolution  of  the  cell  wave 
function  with  each  delta  function  in  the  lattice  simply  places  an  appropriately  phased 
copy  of  the  cell  wave  function  in  every  other  cell  of  the  structure.  The  cell  wave 
function  itself  is  composed  of  a  forward  traveling  or  /  wave  and  a  backward  traveling 
or  g  wave: 

J  0  2  <  -hl2,z>  h/2  .  . 

V’(-)  I  j^jkz  ^  g^-jkz  |,|  <  /j/2. 

Equation  2.1  shows  the  Bloch  wave  function  to  be  a  string  of  identical  but  delayed 
copies  of  the  compound  wave  field  that  exists  in  a  single  cell. 


The  Bloch  wave  function  is  completely  determined  by  q  and  gj /.  While  q  is  found 
from  the  previously  derived  and  experimentally  verified  dispersion  relation  (Eep  2.1 
in  89-6),  g/f,  the  relative  amplitude  (and  phase)  of  the  /  and  g  waves,  is  found  in 
terms  of  impedances: 

/  r  ^Ba  ^oa 

9lf  =  7 — TV~' 

^Ba  +  ^oa 

Here  Zsa  and  Zoa  are  the  acoustic  impedances  of  the  waveguide  with  and  without 
the  side  branches,  respectively.  Note  that  Zbo  depends  upon  frequency  and  various 
waveguide  dimensions.  The  acoustic  pressure  field  measured  at  two  points  in  a  cell 
provided  sufficient  information  to  find  gj /.  Figure  2.1  shows  that  the  theoretical  and 
experimental  results  are  in  excellent  agreement. 


Next  consider  the  nonlinear  case.  Because  the  Bloch  wave  function  is  composed 
of  compound  wave  fields,  our  use  last  year  of  Korpcl’s  progressive  wave  equation  as 
the  model  equation  for  finite  amplitude  Bloch  waves  was  not  really  approi)riatc.  We 
therefore  shifted  this  year  to  the  generalized  Westervclt  c([uation,''  which  is  valid  for 
compound  flow.  A  (luasilinear  analysis  was  employed  in  which  the  fundamental  field 
is  represented  by  its  convolution  expression.  'I'lie  solution  of  the  equation  for  the 
second  harmonic,  found  by  using  a  di.scretc  Green’s  function,  is  as  follows: 

=  G„(.-)  <■  ><(=  -  (i-i) 

n  ——co 

where  qi  is  the  Bloch  wave  number  for  the  fundamental.  The  function  (l„{z)  is  the 
second  harmonic  respon.se  due  to  (he  presence  of  the  fundamental  field  in  a  single  cell. 
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FIGURE  2.1 

The  magnitude  and  phase  of  <^//. 
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i.e.,  the  response  due  to  ^’(2),  and  is  composed  of  both  forward  and  backward  traveling 
Bloch  wave  components.  The  excitation  of  bidirectional  second  harmonic  Bloch  waves 
by  the  fundamental  Bloch  wave  accounts  for  the  large  on-source  second  harmonic  level 
that  was  measured  last  year.  The  present  theory  is  thus  an  improvement  over  the 
progressive  wave  theory. 

A  model  incorporating  the  effects  of  truncation  of  the  periodic  waveguide  and 
second  harmonic  generation  in  the  source  region  was  developed,  and  an  analytical, 
closed  form  result  was  found  (the  series  in  Eq.  2.3  can  be  summed).  The  functional 
form  of  Go{z)  was  derived  by  eissuming  that  the  locally  generated  second  harmonic  due 
to  nonlinear  interaction  of  the  fundamental  /  and  g  waves  is  negligible.  A  comparison 
of  the  theoretical  result  and  a  measurement  is  shown  in  Fig.  2.2.  The  agreement  is 
seen  to  be  fairly  reasonable  (90-3,  90-8). 

Bradley’s  work  seems  to  be  the  first  that  has  ever  been  done  on  nonlinear  acous¬ 
tical  Bloch  waves.  The  only  other  research  of  which  we  are  aware  on  nonlinear  Bloch 
waves  is  in  nonlinear  optics®  and  is  generally  for  the  case  of  small  dispersion.  Disper¬ 
sion  is  definitely  not  small  for  the  Bloch  waves  in  Bradley’s  system. 


2.2  Scattering  of  Sound  by  Sound 


The  purpose  of  this  project,  which  is  3en  Cate’s  Ph.D.  topic,  is  to  perform  an  ex¬ 
perimental  check  of  the  'fjptta’s  predictions  concerning  scattering  of  sound  by  sound. 
This  year  Ten  Cate  concentrated  on  getting  the  new  tank  facility  in  the  Department 


Distance  (m) 

FIGURE  2.2 

Amplitude  of  the  second  harmonic  Bloch  wave  as  a  function  of  distance. 
Fundamental  Bloch  wave  frequency  is  680  Hz  and  source  level  is  133  dB. 
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of  Mechanical  Engineering  (see  Sec.  2.3  in  89- G)  ready  for  nieasureinents.  However, 
together  with  the  Tj0ttas,  he  also  submitted  a  paper,  which  is  a  complete  account  of 
some  earlier  reports  (89-8,  89-9),  to  J.  Acoust.  Soc.  Am.  (90-8). 

The  paper  with  the  Tjottas,  “Effects  of  boundary  conditions  on  the  nonlinear 
interaction  of  sound  beams,”  has  relevance  to  Ten  Cate’s  upcoining  scattering  of  sound 
by  sound  measurements.  In  practically  all  treatments  of  radiation  and  propagation  of 
finite-amplitude  waves,  the  second  order  field  is  presumed  to  be  zero  at  the  sourco(s) 
of  the  primary  field(s).‘‘  The  theoretical  and  numerical  work  in  the  present  paper 
shows  that  the  actual  (nonzero)  boundary  conditions  on  the  second  order  field  have 
effects  that  are  not  negligible.  Figure  2.3  shows  computed  sum  (top)  and  difference 
(bottom)  frequency  fields  for  an  interaction  that  results  in  the  scattering  of  sound 
by  sound.  The  geometry  of  the  transducers  and  the  crossed  beams  they  geneiate  is 
shown  by  the  insert  in  the  center  of  the  figure.  The  boundary  conditions  in  question 
are  those  assumed  on  the  horizontal  dashed  line,  which  [('presents  the  plane  of  the 
sources.  The  dash-dash  dot  patterns  that  peak  at  15°  (top  left  and  lower  right) 
represent  the  classical  product  directivity.  'I'lie  separate  curves  shown  in  the  insets 
show  the  effects  of  different  boundary  conditions  impo.sed  on  the  second  order  field. 
The  four  boundary  conditions  are  as  follows:  (1)  rigid  wall,  (2)  rigid  wall  with  source 
motion  included,  (3)  pressure  reWisc  ;iiul  (1)  pressure  rc'lea.se  \v<tll  with  source 
motion  included.*  'flie  results  show  that  (1)  it  is  mostly  the  scattered  sound  tliat 
is  affected  by  changing  the  boundary  conditions,  and  (2)  the  difference  caused  by 

fifth  curve  is  for  a  nonphysical  bouinlary  coiulilion,  imhnlecl  in  older  to  refi'ieme  pieviou.s 

work.'' 
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FIGURE  2.3 

Beam  patterns  of  scattered  sound  for  various  source  boundary  conditions. 
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Level  (dB  re  Pq  ) 


varying  the  boundary  conditions  seems  to  be,  at  most,  5  to  6  dB.  Also  discussed  in 
the  paper  is  the  relevance  of  boundary  conditions  to  second  harmonic  generation  and 
to  parametric  receiving  and  transmitting  arrays. 

The  rest  of  Ten  Cate’s  time  was  spent  with  the  new  tank  facility  in  the  Mechanical 
Engineering  Department.  The  custom,  computer-controlled  positioning  system  was 
assembled  and  put  in  place  over  the  tank.  The  receive  system  includes  computer- 
controlled  x,  j/,  z  (vertical),  and  B  (about  z  axis)  movement.  The  system  is  closed  loop 
(with  feedback  to  hold  desired  position).  The  transmit  system  includes  a  completely 
manual  control  but  with  computer  readout  of  position. 

Although  the  positioning  system  was  mostly  complete  when  received,  some  work 
was  required  to  make  it  usable.  Two  types  of  receiver  PVDF  hydrophones  will  be 
utilized,  needle  type  and  membrane  type.  Some  machining  had  to  be  done.  The 
membrane  hydrophones  required  a  somewhat  specialized  mount  to  provide  some  tilt 
capability.  The  ultrasonic  sources  required  individually  machined  holders..  The  motor 
drives,  power  sources,  and  associated  equipment  had  to  be  rac':  mounted  and  the 
entire  system  checked  out.  Connecting  and  interfacing  the  Macintosh  II  to  the  motor 
controller,  Berkeley  Axis  Machine  or  BAM,  took  some  time  as  well.  The  BAM  is  a 
computer  itself  and  its  programming  language  had  to  be  learned.  Work  is  currently 
being  done  to  make  the  software  package  LabVIEW  communicate  with  the  BAM  and 
other  instruments  that  will  be  used  in  the  experiments. 

Ten  Cate  also  spent  considerable  time  “tuning”  the  motion  of  each  of  the  axes 
after  everything  was  assembled  and  functioning.  He  aimed  for  repeatability  on  the 
order  of  100  /xm.  Repeatability  is  a  meeisure  of  how  close  one  can  return  to  a  given 
position.  For  reference,  the  PVDF  membrane  hydrophone  used  in  the  tank  has  a 
thickness  of  25  /xm.  Electronic  damping  and  stiffness  are  added  to  each  axis  for 
the  best  motion  profile,  given  the  mass  of  the  sources  and  the  system  parts  being 
moved.  Unfortunately,  the  method  for  arriving  at  the  correct  damping  and  stiffness 
is  essentially  one  of  trial  and  error.  Ideally,  one  desires  smooth  motion  that  produces 
quick  and  accurate  arrival  at  a  chosen  position. 

It  is  expected  that  the  actual  experiments  will  be  well  under  way  before  the 
end  of  the  present  report  period.  Preliminary  qualitative  experiments  are  just  now 
being  started.  Three  types  of  experiments  are  planned:  (1)  reexamination  of  second 
harmonic  generation  from  a  circular  piston  source,  (2)  examination  of  the  effects  of 
pulse  length  on  a  Gaussian  source,  and  (3)  measurement  of  the  scattering  of  sound 
by  sound. 


8 


2.3  Laboratory  Upgrade 


The  first  goal  in  the  upgrade  of  our  Nonlinear  Acoustics  Laboratory  was  to  assem¬ 
ble  a  basic  computer  controllable  data  acquisition  system.  This  task  was  completed 
last  year  (see  Sec.  2.5  in  89-6).  The  heart  of  the  system  is  a  realtime  digitizer  (Tek¬ 
tronix  RTD  710A),  an  arbitrary  function  generator  (Wavetek  model  275),  software  to 
acquire,  process,  and  display  data  (National  Instruments  LabVIEW),  and  a  general 
purpose  interface  bus  (GPIB).  The  last  two  items  enable  our  Macintosh  II  computer 
to  exercise  control  over  the  instruments.  This  year  we  added  important  components 
to  the  system;  a  bandpass  filter  (Krohn  Hite  model  3202),  an  audio  power  amplifier 
(Hafler  Pro  500),  and  a  high  voltage  dc  power  supply  (Classman  EII20P).  Two  more 
items,  a  powerful  spectrum  analyzer  (Il^'wlett  Packard  35C60A)  and  an  A-D/D-A 
board  (National  Instruments  NB-M10-16XII-42),  will  be  purchased  by  the  end  of  the 
report  period.  The  A-D/D-A  board,  to  be  added  to  the  Mac  II,  will  increase  the 
dynamic  range  of  the  entire  system  to  96  dB  and  perform  certain  other  operations. 
Not  all  of  the  instruments  are  used  in  every  experiment,  of  course.  For  example,  the 
high  voltage  supply  is  intended  only  for  our  investigations  with  sparks.  Each  instru¬ 
ment  is,  however,  capable  of  being  an  integral  part  of  any  given  computer-controlled 
measurement. 

We  also  began  the  development  of  a  small  water  tank  facility  for  experiments  in 
the  megahertz  frequency  range.  This  facility  is  not  intended  to  compete  with  the 
larger,  more  elaborate,  and  high  precision  tank  facility  recently  constructed  in  the 
Mechanical  Engineering  Department  on  the  main  campus  of  the  University  (see  Sec. 
2.3  in  89-6).  Since  the  small  tank  will  be  located  at  ARL:UT,  it  will  provide  us  a 
more  convenient  means  with  which  to  carry  out  small-scale  experiments  with  high 
frequency  beams  on  such  topics  as  reflection  and  transmission  at  an  interface.  (It 
will  be  recalled  that  Cotaras’s  study  of  reflection  and  transmission  was  limited  to 
theoretical  calculations,  largely  because  an  airborne  experiment  proved  too  difl[icult; 
see  Sec.  2.1  in  89-6.) 

The  tank  itself  has  been  constructed.  Designed  and  built  of  1/2"  Incite  by  Van 
Doren,  it  has  external  dimensions  23.5"  x  12.5"  x  12.2".  Three  threaded  holes  are 
provided  for  side  or  bottom  mounting  of  transducers  or  spark  sources.  A  spark  source 
has  also  been  constructed,*  and  a  PVDF  needle  hydrophone  has  been  purchased. 

The  proposal  for  a  continuation  of  the  grant  includes  plans  for  equipping  the 
tank  with  hydrophone  positioning  apparatus  and  necessary  electronics  for  sending 
and  receiving  signals.  It  is  hoped  that  the  positioning  apparatus  can  also  be  used 
for  positioning  a  tnicrophone  in  our  airborne  experiments  with  siJarks.  Initially  we 


'Funds  for  building  the  spark  source  came  from  NIil.  'I'lie  tank  may  also  l)c  used  for  .some 
experiments  relevant  to  medical  ultrasonics. 
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had  expected  to  obtain  the  positioning  apparatus  during  the  present  year  (see  Sec. 
2.5  in  89-6),  but  a  one-year  postponement  became  attractive  when  we  learned  that  a 
highly  skilled  person,  Youseph  Yazdi,  will  be  joining  our  group  eis  a  graduate  student 
in  September  1990.  Yazdi  has  had  a  great  deal  of  experience  in  designing  and  using 
tanks  of  this  sort,  and  his  expertise  should  be  invaluable  in  completing  the  facility 
next  year. 


2.4  Miscellaneous 


The  project  on  ellipsoidal  focusing  (see  Sec.  2.4  in  89-6)  continued  at  a  very  slow 
pace  because  of  lack  of  a  student  this  year.  Our  consultant,  Wayne  M.  Wright,  spent 
some  time  rebuilding  the  microphone  preamplifier,  which  had  fallen  on  hard  times, 
and  improving  on  some  of  the  measurements  reported  last  year  (89-2,  89-3).  lie  will 
complete  the  draft  of  a  journal  article  on  the  ellipsoidal  focusing  measurements  before 
the  end  of  the  present  report  period. 

Other  activities  included  presenting  papers  on  evanescent  waves  (90-1),  laboratory 
demonstrations  in  nonlinear  acoustics  (90-2),  absorption  of  finite-amplitude  sound 
(90-6),  and  finite-amplitude  waves  in  porous  material  (90-7). 
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3.  SUMMARY 


During  the  current  report  period,  1  October  1989  -  30  September  1990,  we  have 
been  occupied  primarily  with  two  projects,  (A)  propagation  in  a  periodic  waveguide, 
and  (B)  scattering  of  sound  by  sound.  Work  has  also  been  done  on  (C)  upgrade  of  the 
Nonlinear  Acoustics  Laboratory  and  (D)  some  miscellaneous  items.  Project  (A)  is  a 
combined  theoretical  and  experimental  study  of  acoustical  Bloch  waves  in  a  rectan¬ 
gular,  air-filled  waveguide  that  is  periodically  loaded  with  reactive  branch  elements. 
Although  a  small-signal  study  of  this  system  had  been  completed  the  previous  year, 
discrepancies  between  theory  and  experiment  for  finite-amplitude  waves  prompted  a 
review  of  the  linear  theory.  The  review  led  to  a  new  theoretical  expression  for  the 
nonlinearly  generated  second  harmonic  Bloch  wave.  The  theoretical  prediction  is  now 
in  better  agreement  with  our  measurements.  Project  (B),  although  mainly  experi¬ 
mental,  produced  a  journal  article  this  year  on  the  effect  of  various  source  boundary 
conditions  on  the  scattered  field.  The  Mechanical  Engineering  Department  tank  fa¬ 
cility,  where  the  experiments  on  scattering  will  be  carried  out,  has  been  completed. 
A  series  of  three  different  experiments  is  planned  to  begin  before  the  end  of  the 
present  report  period.  The  main  goal  of  Project  (C),  to  make  possible  a  variety  of 
computer-controlled  experiments,  was  achieved  last  year.  Additional  equipment  was 
purchzised  this  year  to  expand  the  scope  of  the  experiments  that  can  be  done.  Work 
was  also  begun  on  a  small  water  tank  for  experiments  in  the  megahertz  frequency 
range.  Miscellaneous  items  under  Project  (D)  include  work  on  ellipsoidal  focusing 
and  presentation  of  various  papers. 
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